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Abstract 

The interaction between transverse magnetic domain walls (TDWs) in planar 
(2D) and cylindrical (3D) nanowires is examined using micromagnetic simu¬ 
lations. We show that in perfect and surface deformed wires the free TDWs 
behave differently, as the 3D TDWs combine into metastable states with aver¬ 
age lifetimes of 300ns depending on roughness, while the 2D TDWs do not due 
to 2D shape anisotropy. When the 2D and 3D TDWs are pinned at artificial 
constrictions, they behave similarly as they interact mainly through the dipolar 
held. This magnetostatic interaction is well described by the point charge model 
with multipole expansion. In surface deformed wires with artihcial constrictions, 
the interaction becomes more complex as the depinning held decreases and dy¬ 
namical pinning can lead to local resonances. This can strongly inhuence the 
control of TDWs in DW-based devices. 
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1. Introduction 

Magnetic domain walls (DW) dynamics in nanowires is an intensive re¬ 
search t^ic nowadays, with proposed applications ranging from logic to storage 
devices [l|,@,y|. These devices usually exploit the displacement and the pinning 
by artihcially engineered defects of DWs in planar nanowires (nanostrips). 

A large variety of geo metrical constrictions and/or imperfection were studied 
as pinning centers i 1, S 0 to manipulate the DWs position. All types of 
imperfections create a pinning potential for the DW that depends on the DW 
type (transverse or vortex) Q . In the case of a transverse wall, the constrictions 
create a potential well. 
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Transverse walls represent the stable state in low width/diameter nanowires. 
TDW in cylindrical nanowires display different properties than in nanostrips: 
the Walker limit is completely suppressed, theirs velocity and precession speed 
depend linearly on the applied currently while they propagate without deform¬ 
ing theirs internal structure. 

In this paper, we compare the dynamics and the interaction of TDW in ideal 
and rough cylindrical and planar nanowires. In cylindrical ideal nanowires, the 
DWs can form metastable states due to the absence of perpendicular anisotropy 
facilitating the rotation of DWs magnetization in the wall plane[l^. Small sur¬ 
face roughness increases the lifetime of the metastable states due to the local 
pinning. In planar nanowires, the shape anisotropy maintains the DWs mag¬ 
netization in the plane and no metastable state can form. When the DWs are 
pinned at symmetric notches, they interact mainly through the dipolar field in 
both types of wires, and this interaction can be well described by the point 
charge modelfllj. In surface deformed wires, the depinning field decreases and 
DW resonances can occur due to the pinning landscape. The influence of the 
temperature is discussed employing the stochastic ID model. Our results show 
the importance of understanding the pinning mechanism and of mutual inter¬ 
action between DWs in nanowires. 


2. Model and geometry 


We computed a large number of 3D micromagnetic simulations using the 
nmag package [l^ to determine the dynamics of DWs in nanowires 13|. Fig. [I] 
shows the different simulated structures: cylindrical wires with length of 900nm 
or 1200nm and 60nm in diameter and planar wires with same lengths (900nm, 
1200nm) and a width of 60nm and 5nm thickness. The wires were discretized 
into a mesh with a cell size of 4nm, inferior to the exchange length (~5nm 
for Ni). We use the material parameters of Nickel: saturation magnetiza¬ 
tion /roMs=0.6T, exchange stiffness A = 1.05 x 10“^^J/m, 7=188.5GHz/T 


(g factor of 2.15), damping parameter a=0.015 [lj,ll3. No magnetocrystalline 
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Figure 1: (Color online) Simulated structures: (a) perfect cylinder, (b) cylinder with two 
radial constrictions, (c) perfect planar wire and (d) planar wire with two symmetric double 
notches. 
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Figure 2: (Color online) Interaction between three DWs in a perfect Nickel nanowire (length 
= 1200nm, diameter = 60nm): (a) Spatial displacement of DWs during 350ns in a cylindrical 
wire, (b) Synchronization of two DWs from panel (a) determined by the instantaneous period. 

(c) Rotation of the magnetization in the yz plane at the DW center for each DW from panel 
(a). The squares denote the initial position, while the triangles represent the final position. 

(d) Spatial displacement of DWs and (e) precession angle i/; during 400ns in a planar wire. 


anisotropy is considered and the temperature is set to T = OK. The stable struc¬ 
ture in our nanowires is a V-shaped TDW separating domains pointing along 
the long axis (x axis). The average position of each DW was extracted along 
with the azimuthal angle tp of the magnetization inside the DW. 

The constrictions in both types of wires have the same length (20nm) and 
depth (lOnm). In the cylindrical wire they were created as a sharp modulation 
of the nanowires diameter, while in the planar wire the usual symmetric rect¬ 
angular notches are used. The surface deformation was created from a random 
profile along the x axis in both cases, with a predefined amplitude deformation 
a and variable correlation length A. The profile is symmetric on the upper and 
lower edges of the planar wire [3; [I3| and full revolved around x-axis for the 
cylindrical wire. 
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3. Numerical results 


3.1. Perfect wires 

We start with the study of the interaction between DWs in the same ideal 
nanowire in the absence of pinning. The comparison between the two types of 
wires is detailed in Fig. The simulation is initialized with three DWs in the 
nanowire (length 1200nm, diameter/width 60nm), with exterior DWs pointing 
in the same direction and central one pointing at 90° (cylinder) or 180° (strip), 
that are free to relax. In the cylindrical wire (panels (a)-(c)), the magnetization 
in the DW can precess freely around the longitudinal x-axis as no perpendicular 
anisotropy is present. The dynamics in the wire is dominated by the interac¬ 
tion between the DWs. The exterior DWs also interact with the ends of the 
nanowire[l(ij. but this interaction is small as the distance between DW and the 
ends is at least of 340nm, superior to the initial distance between DWs (260nm). 
As initially the DWs are distanced by more than 200nm, the main contribution 
to the interaction comes from dipolar energy which tries to keep their magne¬ 
tization direction antiparallel. The middle DW magnetization rotates to keep 
it antiparallel with the exterior DWs. As they precess, the DWs also oscillate 
spatially along x-axis forming a metastable state (described in detail in 0). 
After 280ns, the left and middle DWs synchronize as can be determined from 
their instantaneous period calculated from the Hilbert transform. The synchro¬ 
nization is a result of the exchange energy which becomes the dominant factor 
in the interaction energy at close range. When the DWs become parallel, they 
can annihilate (320ns) as the exchange energy drops abruptly and the dipolar 
energy remains almost constant. As long as the magnetization in the two DWs 
is antiparallel, they cannot annihilate due to the important exchange energy at 
close separation. This can be observed in panel (d) and (e) for planar nanowires 
with consecutive DWs antiparallel. The DWs cannot annihilate freely as the 
magnetization is frustrated (by the anisotropy) and cannot turn unimpeded in 
the yz plane. Initially, the DWs attract due to the dipolar field and they ap¬ 
proach each other, but rapidly after 2ns they repel and remain stationary at 
a fixed distance. The ip angle (panel (e)) stays along the y direction having a 
small initial oscillation (^ 3°). If the DWs are initiated with the same direction 
of magnetization in planar nanowires, they annihilate rapidly in less than Ins 
(image not shown) [l8| . 

3.2. Rough wires 

Usually, nanowires always present some kind of surface roughness. The 
cylindrical wires were radially deformed while the planar wires were edge de¬ 
formed. Two deformation amplitudes were used: cr=2nm with A=3-10nm and 
(T=4nm with A=5-15nm. An example of meshes used in simulations are shown 
in Fig. (Ha). 

To compare with the ideal case we used the same initial configuration of three 
DWs situated at equal distance. Five different meshes were studied for each 
deformation amplitude. We observe two different behaviors in cylindrical wires: 
at small deformation (as in Fig. [3Kb) and (c)), the DWs sometimes effectuate 
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Figure 3: (Color online) Meshes of the surface deformed wires used in simulations (a). Spatial 
displacement of DWs and theirs instantaneous period for a cylindrical nanowire with radial de¬ 
formed surface as in (a): (b) and (c) correspond to soft roughness with amplitude deformation 
(T=2nm (upper wire in (a)) while (d) and (e) to hard roughness ((j=4nm). 


small jumps depending on the local pinning profile and annihilate after a longer 
time than in the smooth case (around 400ns) for four cases out of five. In the 
case presented in panels (b) and (c), the left and middle DW rotate clockwise 
initially. After 100ns, the middle DW starts to change its rotation direction 
and to rotate in the same direction as the right DW (anticlockwise rotation). 
It oscillates around a fixed position with a small amplitude (8nm) for almost 
200ns, while the other two DWs slowly move closer. The right DW moves 
closer to the middle DW (which stays pinned) compared to the left DW and 
therefore it interacts stronger (the angular velocity of the DW is proportional 
with the magnetostatic fieldQ). The magnetization in the middle DW starts 
to oscillate between clockwise and anticlockwise direction until finally stabilizes 
and starts to synchronize with the right DW. When it finally synchronizes, after 
380ns, it moves away form its local trapping position and annihilates. At higher 
deformation (Fig.jSjd) and (e)), the DWs sometimes effectuate larger jumps and 
annihilate after an identical or lower time than the ideal case (around 300ns) for 
four cases of five trials. In the presented case, the left and middle DWs rotate 
clockwise, but rapidly after 50ns, the middle DW starts to change its rotation 
direction as it jumps closer to the right DW while the lower DW remains pinned 
at a local position. After 150ns, it oscillates almost synchronously with the 
right DW but it does not annihilate until 320ns due to the larger spatial jumps. 
Local pinning sites trap the DWs and a larger interaction energy is needed to 
depin them hindering the annihilation. This local pinning also influences the 
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Figure 4: (Color online) (a) Equilibrium position for a pinned DW function of the applied dc 
magnetic field for a planar or cylindrical wire with a symmetric constriction at its center. The 
oscillations of the spatial position and azimuthal angle of the DW are shown for a dc field of 
290Oe for a cylindrical wire in (b) and for a dc field of 490e for a nanostrip in (c). 


interaction between DW pairs, some DWs can approach while others stay pinned 
and therefore can change theirs rotation direction. 

The edge roughness does not influence the interaction of DWs in planar 
nanowires, as the DWs magnetization is still trapped in-plane. Therefore, no 
metastable states appear, with or without surface deformation. 

3.3. Pinning at artificial constrictions in a perfect wire 

To be able to compare the DWs interaction between the two types of wires, 
we examined the case when the DWs are pinned at artificial constrictions in 
a smooth nanowire. For the planar wires we used two symmetric notches as a 
pinning center, while for the cylindrical wires we used a radial constriction (as 
shown in Fig. [T|). 

We first evaluate the static interaction of a DW with a single pinning center, 
situated in the center of the wire, from quasistatic micromagnetic simulations. 
The average equilibrium position is shown in Fig.lJJa) for both wire types. For 
planar wires two symmetric rectangular notches create a harmonic potential 
(x linear in JI), while a radial notch in a cylinder creates a complex potential 
with a change of slope when the DW center is pinned outside the notch[l9l|. 
The calculated spring constants and depinning fields are k = 8x 10“®N/m and 
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Hdep=50Oe for the planar wire, fci = 1.23 x 10“^N/m and /c 2 = 3.21 x lO^'^N/m 
and Hdep=3000e for the cylindrical one. The radial notch has a larger depinning 
field and the DW relaxes to equilibrium doing several rotations in the yz plane 
(Fig. SKb)), however the DW pinned at a radial notch has also a very long 
relaxation time (hundreds of ns) compared to the symmetric rectangular notches 

(Fig. He)). 

The domain wall width suffers a small modification when submitted to the 
external field. The width contracts by 2 nm just before depinning from 20nm 
(without field) to 18nm for the cylindrical wire and from 31.9nm to 30.9nm for 
the nanostrip. 

We analyze next the DWs interaction considering two DWs of opposite chi¬ 
rality head-to-head and tail-to-tail (HHDW and TTDW) pinned at notches 
situated in the same wire. We start with the cylindrical wire, where several 
edge-to-edge distances between notches were used: 180nm, 230nm and 270nm, 
with identical notch form. The DWs interact through the dipolar field and 
theirs equilibrium positions in absence of external perturbation is not at the 
center of each notch. They precess around the cylinder axis and relax with the 
magnetization pointing in opposite directions at theirs center, corresponding to 
an attractive interaction, to minimize energy. An example of the interaction 
between the two DWs is shown in Fig. [Sj8')"(b) for the 180nm separation dis¬ 
tance, where the orientation of the DWs center magnetization in the yz plane is 
shown. An external magnetic field of -l-lOOe is applied to push closer the DWs. 
Theirs initial orientation, indicated by squares, is opposite and after 200ns they 
relax to equilibrium maintaining an opposite orientation of magnetization (final 
position indicated by triangles). In panels (c) and (d), the relaxation process 
is shown for the right DW and a very small oscillation around the equilibrium 
position is visible after the 200ns due to the local pinning potential as the center 
of the DW is situated at the edge of the notch. The closest position of the two 
DW before depinning is illustrated in the inset of the panel (f) for an applied 
field of -l-1250e. 

Due to their interaction, the equilibrium position of each DW is closer to the 
other notch when the notches are close together, as can be observed in Fig. [SJe) 
where is represented the potential profile variation with the distance between the 
two DWs. The minimum of the potential well (no external perturbation present) 
is outside of the notches for the closest separation (180nm), when the DWs are 
situated at a distance of 13.5nm from the center of each notch. The pinning 
potential is highly distorted in this case as more field is needed to separate them. 
The interaction between the DWs also greatly changes the depinning field for 
the two directions of the nanowire. 

To quantify the interaction between the two DWs pinned at notches, we 
represent in Fig. [SJf) the variation of the interaction energy between the DWs 
with their separating distance. The interaction energy was calculated as follows: 
first the energy of the nanowire without any DW was calculated (edge energy). 
Then, the pinning energy for each DW was calculated similarly with the case 
of a single notch: we initialized the nanowire with a single pinned DW (left or 
right) and determined the potential profile of each of the notches. Afterwards, 
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Figure 5: (Color online) Interaction between two DWs pinned at notches: (a) and (b) Rotation 
of the magnetization in the yz plane at the DW center for two DWs pinned at notches separated 
by 180nm, when an external field of lOOe was applied. The squares denote the initial position, 
while the triangles represent the final position, (c) and (d) Relaxation of the DW position 
and precession angle for the DW from (b). (e) Variation of the normalized potential energy 
of the pinned DWs with the separation distance between them for three notch separation of 
180, 230 and 280nm. (f) Variation of the interaction energy between the two DW with the 
separation distance (circles) for the same three notch separation as in (e). The line represents 
a fit with the multipole expansion. The inset shows two DWs pinned at notches separated by 
180nm. 


we extracted from the total energy the pinning energy of each DW and the initial 
demagnetizing energy (edge energy), to obtain the interaction energy. In the 
figure, we represented only the interaction energy (normalized) when pushing 
the DWs closer for each of the three notch separations studied. As observed 
from the fitted curve, the calculated points fall on an universal curve which 
correspond to the dipolar interaction (multipole expansion) as discussed below. 

The interaction energy between two DWs of opposite chirality in the same 
wire was calculated using the point-charge model with multipole expansion 111. 
The expression disregards the internal structure of the DWs and the exchange 
interaction between them (thus applicable at high separation of DWs), consider¬ 
ing only the interaction through the dipolar field therefore giving a qualitative 
description of the interaction. The interaction energy between two DWs of 
opposite chirality (HHDW and TTDW) separated by d is: 


E = 


Mogi'?2 


47rd 
Trd 


1 - 


tt^A^ cos('(/'i - V'2)) 


1 - 


4d2 

tt^A^ cos(-!/'i - '02)) 
4d2 


^Qq 


-E 


I2d? 


( 1 ) 


where qij represent the ’magnetic charges’ of the DWs and in this approxi¬ 
mation are equal to ±2MsS, with S the section of the nanowire and A the 
DW width. The first term represents the Coulomb-type interaction energy 
(monopole-monopole), while the second is the correction from the dipole-dipole 
interaction. The third term represents the quadrupole-monopole contribution 
with w the width and t the thickness of the nanowire. We fitted the above 
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Figure 6: (Color online) (a) Variation of the normalized potential energy of the pinned DWs 
in nanostrips with the separation distance between them for three notch separation of 180, 
230 and 280nm. Two cases for each separation are shown: DWs magnetization parallel or 
antiparallel, (b) Plot of the interaction energy between two pinned antiparallel DWs. 


expression to the numerical values of Fig. [5)f) (line). From the parameters of 
the fitted curve, we extracted the domain wall width and obtained a value of 
23.28nm, larger than the value obtained by fitting the profiles of DWs with the 
Bloch profile which varies between 17.91nm and 19.67nm depending on the ap¬ 
plied field and local pinning position. This discrepancy between the two models 
comes mainly from the error in estimating the average position of the DW in the 
3D wire as the DW distorts taking a bow shape when pinned at the notch while 
submitted to an external field. This general dependence of the interaction of 
DWs describes the importance of the dipolar energy when the DWs are situated 
at large distances, where the exchange energy is much smaller (here by a factor 
of five). 

The same analysis as above was effectuated for two DWs of opposite chirality 
pinned at symmetric rectangular notches inside the same nanostrip (Fig. [TJd)). 
The DW interaction was studied for several distances between notches (180nm, 
230nm and 280nm). As the magnetization of the two DWs is confined to the 
y direction, two cases were studied: parallel magnetization of the DWs or an¬ 
tiparallel. The potential energy of the pinned DWs varies with the separation 
distance and is shown in Fig. [HKa) for three notch separations. 

When the DWs are pinned at larger distances (230nm, 280nm), the potential 
energy is quasi identical for parallel or antiparallel DWs. The most important 
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Figure 7: (Color online) Variation of the equilibrium position, with applied field, of two pinned 
DWs for ISOnm notch separation in a cylindrical wire (a)-(b) and respectively in a planar wire 
(e)-(f). The ideal case is compared with several surface deformed cases as explained in the 
text. The variation of the equilibrium position for 230nm notch separation is shown in (c)-(d) 
for the cylindrical wire and in (g)-(h) for the planar wire. 


difference appear at ISOnm notch separation. When the DWs are parallel, a 
smaller positive depinning field (DWs pushed closer) is needed as they attract 
each other and after depinning will annihilate. When the DWs are antiparallel, 
the positive depinnig field is larger as the DWs repel each other due to the large 
exchange energy at close separation. After depinning the DWs don’t annihilate 
but form a 360° DW[l8|. The anisotropy impedes the magnetization to rotate 
inside the DWs and the antiparallel DWs cannot annihilate. The interaction 
energy between DWs was extracted in the same manner as for cylinders and is 
shown in Fig. |6l(b). A DW width of 29.30nm was determined from the point 
charge model, that is closer to the calculated Bloch width (30.96nm - 31.95nm) 
than in the 3D case. In 2D wires, the error when extracting the DW average 
position is lower allowing for a better agreement between the calculated DW 
width by the different models. 

3.4- Pinning at artificial constrictions in a rough wire 

The roughness modihes the interaction between DWs pinned at notches. 
Several simulations were computed for surface deformed wires with cr=4nm and 
A=3-15nm, using radial deformation for cylinders (5 cases) or edge deformation 
(10 cases) for strips. The results are summarized in Fig.[3 The variation of the 
equilibrium position of pinned DWs with applied field for ISOnm and 230nm 
notch separation is depicted for both wires in the ideal and rough cases. 

For the deformed cylindrical wire with DWs pinned at ISOnm notch sepa¬ 
ration ((a) and (b)), at positive applied field (DWs pushed closer), the dipolar 
interaction is strongest and prevails the local pinning due to roughness. The 
position of the DWs follows almost the same curve as in the ideal case outside 
of the notch (same slope). For negative applied field, the DW travels through 
the notch and small variations to the ideal case can appear inside the notch 
and at large fields. At 230nm notch separation ((c) and (d)), the influence of 
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the surface roughness is more pronounced and dynamical pinning is more eas¬ 
ily observable. The depinning field decreases in all deformed wire simulations 
compared to the ideal wire, when the DWs are pushed closer (positive field), 
by 25-350e at ISOnm notch separation and by 25-50Oe at 230notch separa¬ 
tion. It also decreases in four out of five cases for negative applied fields with 
same amounts as for positive fields. In one case, for each notch separation, 
the depinning field was identical as the ideal one. In terms of the variation of 
normalized potential energy with the separation distance between DWs (as in 
Fig. [S{e)), the DWs depin more rapidly when pushed together and the height 
of the pinning potential well by the notch is diminished. The potential well 
remains asymmetric with a larger deviation from the ideal case when the DWs 
are pinned inside each notch. The interaction energy between DWs remains 
magnetostatic, but a quantitative determination is more complicated due to the 
lack of precise determination of the pinning energy due to roughness. However, 
the opposite analysis can be conducted and the pinning energy due to roughness 
can be determined quantitatively considering that the interaction energy and 
the pinning energy by a notch are completely determined from the ideal case. 

In deformed nanostrips, a change of slope is observed in eight out of ten 
simulations, compared to the ideal nanostrip. The average position of the DWs 
of all ten rough nanostrips is shown, along with all individual cases, in panels 
(e) to (h) for both notch separation. The average curve (dotted line) has the 
same slope as eight of the rough nanostrips. The two remaining simulations 
follow the ideal case or have a different slope. At ISOnm notch separation the 
slope variation is lower than for the 230nm separation, as in cylinders. The 
depinning field decreases in all cases for 230nm case and for negative fields for 
ISOnm case with 8-140e. For the positive applied fields at ISOnm separation, 
the depinning field decreases in six out of ten cases and in the remaining four 
increases slightly with maximum 60e. 

In the case presented in Fig. IHa) for a deformed cylindrical wire, the right 
DW is pinned at the same equilibrium position situated 15nm from the center 
of the right notch (at 125nm) when the applied field is varied between 50Oe and 
1250e. As can be observed in panels (a) and (b), around 750e there is a reso¬ 
nance of the right DW in the local potential, as the DW position oscillates with 
high amplitude and the magnetization processes in a steady state around the x 
axis. At lower fields, the magnetization only processes initially and afterwards 
the oscillation damps to the equilibrium position. At resonance, a steady state 
develops similar to an applied polarized current when the spin torque counter 
balances the damping. The DW resonance depends on the local potential curva¬ 
ture and therefore on the geometry of the deformation (the deformation profile is 
close to triangular). This dc magnetic induced resonance resembles the current 
induced resonance of geometrically confined DWs in rings and in triangular 
notches in stripes[2l|, although without the transformation in the DW struc¬ 
ture. The DW stays pinned until 1250e, although the spatial oscillation of the 
DW position diminishes. If the field is increased above the local depinning field 
(150Oe), the DW makes a jump of 10 to 15nm and behave like in the perfect wire 
relaxing to a local equilibrium. The resonance phenomena in rough cylindrical 
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Figure 8: (Color online) (a) Plot of the right DW position and (b) azimuthal angle in a 
deformed cylindrical wire, for different applied fields, showing local resonance, (c) Plot of the 
DW position and (d) azimuthal angle in a deformed planar wire for H = 420e. 


wires is observed in all the computations for one or both pinned DWs. 

No resonance is observed in the deformed nanostrips, only an asymmetry in 
the oscillation when relaxing to the equilibrium position as in panels (c) and 
(d) of Fig.[5l 


4. Discussion and Conclusion 


The interaction energy between DWs in the same nanowire was shown to 
be mainly dipolar and depend on the roughness. However, the DWs interaction 
depends also on temperature as the pinning of a DW by a constriction in a 
nanowire can be greatly influenced by it. This effect is often estimated using 
the stochastic one dimensional model which assumes the DW as rigid and is 
usually determined from 22, 2^ : 
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with X and ip the position and respectively the azimuthal angle of the DW, Hk 
the DW demagnetizing field, qx and q^ represent stochastic Gaussian noise with 
zero mean value and correlations {qi{t)qj[t')) = {2akBT)/{fioMsAS)6ijd{t — t'). 
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E is the potential energy of the DW that includes the internal energy, the 
Zeeman energy, the interaction energy with other DW and the pinning energy. 

We evaluated the depinning probability of a DW from a pinning potential 
corresponding to the single radial notch in the cylindrical nanowire or to the 
symmetric rectangular notches in nanostrips. The same parameters as in the 
micromagnetic simulation were used and 1000 realizations were computed for 
each applied field. In the nanocylinders, at T=0K, the depinning occurs at 
2990e, while at T=300K the depinning probability changes from 0% to 100% 
in a field range of 140e (2830e to 2970e). In planar wires, at T=0K, the 
depinning takes place at 50Oe, while at T=300K the depinning is of 3% at an 
applied field of 250e. Based on this analysis, the thermal fluctuations shoud 
impact more the pinning of a DW in nanostrips. This can be deduced also 
from the height of the pinning potential, which is less than O.SeV for notches in 
nanostrips compared with several eV for the radial notch in nanocylinders. 

The cylindrical nanowires have advantages but also drawbacks for practical 
applications as compared to nanostrips. Cylindrical constrictions give much 
longer DW relaxation time, which increases with roughness. However, as no 
Walker limit is present, the study of DW interactions in such structures could be 
the basis for in-depth exploring of fundamental interactions. We speculate that 
triangular constrictions will lead to interesting resonance states in an applied 
dc magnetic field as seen in the rough wires. 

In conclusion, we have presented a systematic study of the interaction be¬ 
tween TDWs in planar and cylindrical nanowires. Free DWs form metastable 
states in free anisotropy wires (as in the cylindrical geometry). Two pinned DWs 
interact strongly through the dipolar field in both wire types. Surface rough¬ 
ness modifies the lifetime of metastable states and the interaction of pinned 
DWs becomes more complex due to the asymmetry, dynamical pinning and 
temperature. 

The authors wish to thank L. Raymond and is grateful for the support of 
the NANOMAG platform by FEDER and Ville de Marseille. The computations 
were performed at the Mesocentre d’Aix-Marseille University. 


References 

[1] D. A. Allwood, G. Xiong, C. C. Faulkner, D. Atkinson, D. Petit, R. P. 
Cowburn, Science 309, 1688 (2005). 

[2] S. S. P. Parkin, M. Hayashi, L. Thomas, Science 320, 190 (2008). 

[3] M. Hayashi, L. Thomas, R. Moriya, G. Rettner, S. S. P. Parkin, Science 
320, 209 (2008). 

[4] M. Klaui, J. Phys.: Condens. Matter 20, 3I300I (2008). 

[5] S. Glathe and R. Mattheis, Phys. Rev. B 85, 024405 (2012). 

[6] D. Petit, A.-V. Jausovec, D. Read, and R. P. Cowburn, J. Appl. Phys. 103, 
II4307 (2008). 


13 



[7] M. Franchin, A. Knittel, M. Albert, D. S. Chernyshenko, T. Fischbacher, 
A. Prabhakar, and H. Fangohr, Phys. Rev. B 84, 094409 (2011). 

[8] D. Petit, A.-V. Jausovec, H. T. Zeng, E. Lewis, L. O’Brien, D. Read, and 
R. P. Cowburn, Phys. Rev. B 79, 214405 (2009). 

[9] M. Yan, A. Kakay, S. Gliga, and R. Hertel, Phys. Rev. Lett. 104, 057201 

( 2010 ). 

[10] V. O. Dolocan, Eur. Phys. J. B 87, 188 (2014). 

[11] B. Kruger, J. Phys.: Condens. Matter 24, 024209 (2012). 

[12] T. Fischbacher, M. Franchin, G. Bordignon and H. Fangohr, IEEE Trans. 
Mag. 43, 2896 (2007) (http://nmag.soton.ac.uk). 

[13] V. O. Dolocan, Appl. Phys. Lett. 101, 072409 (2012). 

[14] K. Gilmore, Y. U. Idzerda, and M. D. Stiles, Phys. Rev. Lett. 99, 027204 
(2007). 

[15] B. Heinrich, Ultrathin Magnetic Structures III (Springer, New York, 2005). 

[16] Y. Nakatani, A. Thiaville, and J. Miltat, Nat. Mater. 2, 521 (2003). 

[17] E. Martinez, J. Phys.: Gondens. Matter 24, 024206 (2012). 

[18] A. Kunz, Appl. Phys. Lett. 94, 132502 (2009). 

[19] V. O. Dolocan, Appl. Phys. Lett. 105, 162401 (2014). 

[20] D. Bedau, M. Klaui, M. T. Hua, S. Krzyk, U. Rudiger, G. Faini, and L. 
Vila, Phys. Rev. Lett. 101, 256602 (2008). 

[21] A. Bisig, M. Klaui, M. T. Hua, S. Krzyk and U. Rudiger, App. Phys. Lett. 
95, 162504 (2009). 

[22] M. E. Lucassen, H. J. van Driel, C. Morais Smith, and R. A. Duine, Phys. 
Rev. B 79, 224411 (2009). 

[23] E. Martinez, L. Lopez-Diaz, L. Torres, C. Tristan and O. Alejos, Phys. 
Rev. B 75, 174409 (2007). 


14 


